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Abstract 

An introduction to dynamical microscopic models of hadronic and nuclear 
interactions is presented. Special emphasis is put in the relation between 
multiparticle production and total cross-section contributions. In heavy ion 
collisions, some observables, considered as signals of the production of a Quark 
Gluon Plasma (QGP), are studied. It is shown that they can only be described 
if final state interactions are introduced. It is argued that the cross-sections 
required are too small to drive the system to thermal equilibrium within the 
duration time of the final state interaction. 
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1 Introduction 



Statistical QCD predicts the existence of new states of matter at high temperature T 
and high baryon number densities fi (baryochemical potential). The phase diagram 
is schematically represented in Fig. 1. Let us discuss first the phase transition at 
high T and small fi (i.e. when baryon and antibaryon densities are approximately 
equal). At n — 0, lattice calculations PQ show a phase transition to a deconfined 
plasma of quarks and gluons (QGP) at a critical temperature T c ~ 150 4- 200 MeV, 
corresponding to an energy density e - an order of magnitude higher than that 
of ordinary nuclear matter (eq = 170 MeV/fm 3 ). The results are shown in Fig. 2. 
Below the critical temperature T c , we heat the system and the temperature increases 
towards T c . At T ~ T c , the temperature remains constant (the energy given to the 
system is used to increase its latent heat) and the energy density e increases sharply. 
For T > T c the temperature increases again with e/T 4 approximately constant, as 
it should be for an ideal gas (Stefan-Boltzmann limit). The sharp increase of e/T 4 
near T c is due to the increase in the number of degrees of freedom of the system 
from a hadronic phase (T < T c ) to a plasma of quarks and gluons (T > T c ) . 

In the ideal case of pure gauge QCD (the so-called "quenched" approximation 
in which dynamical quarks are absent) the phase transition is first order. In the 
presence of dynamical quarks the situation is more complicated and the order of the 
phase transition depends on the number of flavors. A restoration of chiral symmetry 
also takes place at the same T c . 

The region \i ~ studied in lattice QCD corresponds to the conditions in high- 
energy heavy ion collisions at mid-rapidities. It also corresponds to the conditions in 
the primordial universe. In the opposite conditions, i.e. low temperature and high 
baryon densities (i.e. the conditions at the center of neutron stars) a new phase (or 
phases) is expected j2] producing a color superconductor in which pairs of quarks 
condensate - in a way similar to electron (Cooper) pairs in QED. 

A possibility to create in the laboratory the high energy densities and tempera- 
tures required for the production of QGP is via head-on heavy ion collisions at high 
energy. In AA collisions, the energy density can be evaluated using the Bjorken 
formula 



Here dN/dy is the number of produced hadrons per unit rapidity and < my > is 
their average energy. Ra is the nuclear radius and r an average formation time. 
Taking the customary value Tq ~ 1 fm, we obtain for central Au Au collisions at 
RHIC an energy density e ~ 4 GeV/fm 3 , which appears to be high enough for QGP 
production. 




2 



Does it mean that QGP production in high-energy heavy ion collisions follows 
from QCD ? The answer is negative. Indeed, QGP formation is predicted by sta- 
tistical QCD, i.e. QCD applied to a system in thermal equilibrium. Therefore, one 
of the main issues in heavy ion physics is to determine whether the produced final 
state reaches thermal equilibrium. This depends, of course, on the strength and 
time duration of the final state interaction and can only be decided with the help of 
experiment. The observables which provide the most reliable information are the so- 
called signals of QGP. In this lecture I will discuss two of the most important signals : 
the particle abundance (in particular, hyperon and antihyperon production) and the 
J/ip suppression. The latter could be due to the deconfinement or "melting" of the 
cc bound state in the plasma [3] . However, a similar phenomenon is observed in pA 
collisions, where the produced densities are too small for QGP formation. As for the 
former, an argument in favor of equilibrium is the fact that particle abundances are 
well described using statistical models Moreover, these models provide a natural 
explanation [5 of the increase of the relative yields of strange particles in central 
nucleus- nucleus as compared to pp collisions (strangeness enhancement). However, 
one should take into account that these models are also very successful jH] in pp and 
even in e + e~ interactions, where QGP is not produced. 

Therefore, it is important to study these observables in the framework of micro- 
scopic models which are successful in describing pp and pA interactions and can be 
generalized to heavy ion collisions. Models of this type [7-10] are called string mod- 
els, in which particle production takes place in the form of strings (chromo-electric 
flux tubes) stretched between constituents of complementary color charge. In QCD, 
the force between complementary color charges (such as a quark and an antiquark 
of same color), are small at distances smaller than the hadron radius. However, at 
larger distances, the potential of the chromo-electric field increases linearly with the 
distance (confinement). Due to the strong force generated by this potential it is not 
possible to isolate the two color charges. Indeed, in the attempt to separate them, 
the potential energy of the system increases very rapidly and is converted into mass 
via the creation of quark-antiquark pairs. This results in the production of hadrons 
(mostly mesons), formed by recombination of a quark and an antiquark of adjacent 
pairs. The lines of force of the chromo-electric field are strongly colimated along the 
axis determined by the two color charges. Hence the string-like or jet-like shape of 
the set of produced hadrons. 

As a starting point, one usualy assumes in these models that particles produced 
in different strings are independent. In this case thermal equilibrium cannot be 
reached, no matter how large the energy density is. Indeed, in this large 
energy-density is the result of piling up a large number of independent strings. 
In other words, some "cross-talk" between different strings is needed in order to 
thermalize the system. 

The assumption of independence of strings works remarkably well in hh and 
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hA interactions [7-10] - even in the case of event samples with 5 or 6 times the 
average multiplicity - indicating that no sizable final state interaction is present 
in these reactions. However, it is clear that in heavy ion collisions, where several 
strings occupy a transverse area of 1 fm 2 , the assumption of string independence 
has to break down. This is indeed the case. As we shall see, some data cannot 
be described without final state interaction. It could have happened that this final 
state interaction is so strong that the string picture breaks down and becomes totally 
useless. This does not seem to be the case. On the contrary, present data can be 
described using the particle densities computed in the model as initial conditions in 
the gain and loss (transport) equations governing the final state interaction. The 
interaction cross-section turns out to be small (a few tenths of a mb). Due to this 
smallness and to the limited interaction time available, final state interaction has 
an important effect only on rare processes, in particular 5, Q and J/ip production, 
or particle yields at large pr- The bulk of the final state is not affected. 

The plan of these lectures is as follows : in Section EJ I introduce the general 
framework of high-energy scattering in a hadronic language with special emphasis on 
the unitarity condition and its implementation in eikonal and Glauber Models. This 
section contains also some rudiments of Regge poles and the concept of Pomeron. 
In Section |2J I introduce a microscopic model : the Dual Parton Model (DPM) and 
compute the charged particle multiplicities as a function of energy and centrality 
(the latter characterizes the impact parameter of the collision in a way which is 
experimentally measurable). In Section HI I study the so-called stopping power, i.e. 
the fate of the nucleons of the colliding nuclei. Particle abundances are studied in 
Section and J / ip suppression in Section EJ Section [7| contains the conclusions. 

2 High-Energy Scattering : Hadronic Picture 

2.1 General Framework 

An important property of strong interactions is the unitarity of the S matrix oper- 
ator 

ss + = S + S = 1 . (1) 

The S matrix is the operator that transforms free states at time t = — oo into free 
states at times t = +oo. If we write 

S = 1 + iT , (2) 

where < a\T\b > is the transition amplitude between free states a and b, the unitarity 
condition ((TJ) reads 
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2 Im <a\T\b>= ]T < a\T\n > < n\T\b >* . (3) 
all n 

Here n denotes any state that can be reached from both a and b. For a identical to 
b, eq. reduces to 

2 Im <a\T\a>= J2 (< a\T\n >) 2 . (4) 
all n 

A two body amplitude 1 + 2 — > 3 + 4 depends on two independent variables : 

s = (pi+ P2) 2 , t=(pi- p 3 f , u = (pi - p 4 ) 2 , (5) 

with s + 1 + m = + m\ + m| + m^. When a and 6 are identical, t = 0. Since 
the sum in the r.h.s. of (@J) builds up the total 1 + 2 cross-sections, we obtain from 
(0J) the optical theorem, which relates the total cross-section to the forward (t = 0) 
imaginary part of the elastic amplitude T e [(s,t) 

Ait 

crtot(s) = — Im T d (s,t = 0) . (6) 
s 

Let us now examine some of the implications of the unitarity condition It 
is convenient to separate the states n into two classes. The first class consists of all 
inelastic states containing no large rapidity gap between the produced particles. The 
second class consists of those events with at least one large rapidity gap. Obviously, 
elastic scattering produces a state of the second class. We know experimentally that 
the first class of events gives the most important contribution to the total cross- 
section. It is convenient to depict such a contribution in the form of a diagram as 
shown in Fig. 3a. (Its interpretation will be discussed below). The important point 
is that the final state in Fig. 3a has also to be included as intermediate state in 
4. This is precisely the contribution of the elastic state which, as discussed above, 
belongs to the second class of intermediate states. One obtains in this way the 
diagram shown in Fig. 3b. Here the lines with a cross denote the contribution of the 
on-shell initial state to the elastic amplitude. Repeating such an iteration, one is led 
to the multiple scattering diagram depicted in Fig. 3c. One has to add up all the 
diagrams in Figs. 3. In fact, inelastic states (for instance pN* and N*N* in the case 
of pp scattering) also contribute as intermediate states, so that the vertex function 
(blob) of the diagrams in Fig. 3b and 3c can be quite complicated. There are also 
contributions due to intermediate states n containing more than one rapidity gap. 
However, their contributions are small at present energies and will not be considered 
here. 

The interpretation of the diagram in Fig. 3a originates from the claim that the 
inelastic intermediate states (without rapidity gap) generate a Regge pole called the 



5 



Pomeron (see Section I23j) . All diagrams can be treated as Feynman diagrams in a 
field theory called Reggeon Field Theory or Gribov's Reggeon Calculus [llj. When 
only the diagrams in Figs. 3 are kept, one obtains, as a particular case, the so-called 
eikonal model for hadron-hadron scatttering or the Glauber model in interactions 
involving nuclei. 

2.2 The Eikonal and Glauber Models 

In the case of hadron-hadron interactions the diagram in Fig. 3a (Pomeron) is also 
called the Born term. Its contribution to the total cross-section is parametrized 
in the form g-&{t) gzdf) s A where g are coupling constants and A a parameter (see 
Section l23j) . Let us consider the case of pp scattering and let us assume for simplicity 
an exponential dependence, in t, i.e. 

013 (*) =92i(t) = Aexp(Bt) . (7) 

In this case, all loop integrals in Figs. 3b and 3c can be performed analytically. At 
a given s there are only two free parameters A 2 s A and B, which can be determined 
from the experimental values of a tot and a e i (or from a tot and the slope of the elastic 
amplitude). For s variable there is a third parameter A, and the values of <j to t an d 
a e i at several energies can be described in this way. 

A contribution to the total cross-section in s A with A > violates the Froissard 
bound [12] (consequence of unitarity). Note that A > is needed in order to describe 
the rise with energy of a tot . However, when all multiple scattering terms are added, 
one obtains, at high energy, a to t ^ £n 2 s, i.e. the maximal increase with s allowed 
by the Froissard bound. 

It is convenient to work in the impact parameter (b) representation. The scat- 
tering amplitude in 6-space is defined as T(s,b) = (1/2%) J d 2 qr exp(—iqr • b)T(s,t) 
where t = —q\- Using this transformation it is easy to see that at fixed b the con- 
tribution of a diagram involving n exchanges (Fig. 3c) is just the n-th power of the 
Born term (Fig. 3a) times some trivial combinatorial factors. Furthermore the sum 
over n has a simple expression (see below, footnote 1). All these features will now 
be described in detail in the case of the Glauber model. 

Let us consider, for definiteness, a proton-nucleus (pA) scattering. In this case 
there are no free parameters involved. Indeed, the t-dependence (or 6-dependence) 
of the proton vertex function can be neglected in comparison with the fast variation 
of the nuclear one. The latter is known from the nuclear density. Moreover, in this 
case the Born term at t = is just the inelastic proton-nucleon cross-section, which 
is known experimentally. The main formula of the probabilistic Glauber model is 
the one that gives the cross-section a n for n inelastic collisions of the projectile 
with n nucleons of the target nucleus (the remaining A — n nucleon which do not 
participate in the interaction are called spectators), at fixed impact parameter b : 
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a n {b) = Q (a mel T A (b)) n (1 - a mel T A (b)) A ~ n . (8) 

Here Omei is the proton-nucleon inelastic cross-section and TA(b) is the nuclear pro- 
file function (obtained by integrating the nuclear density : 7a (6) = f_ ™ dZp A (Z, b) ; 
/ d 2 bT A (b) = 1). Equation (JHJ) is just the Bernoulli's formula for composite probabil- 
ities. The first factor is a trivial combinatorial factor corresponding to the different 
ways of choosing n nucleons out of A. The second one gives the probability of hav- 
ing n inelastic pN collisions at given b. The third one is the probability that the 
remaining A — n nucleons do not interact inelastically. Let us consider first a term 
with two collisions both of which are inelastic. The corresponding cross-section is 
a 2^P) = {^)( a ineiTA(b)) 2 i.e. a positive term. Let us now consider the case of two 
collisions only one of which is inelastic. The corresponding (interference) term is 
cr\(b) obtained from Eq. (JHJ) by putting n = 1 and taking the second term in the 
expansion of the last factor. We get tj\{b) = —A (A — \){a ine iT A {b)) 2 . We see that 
<y\(b) = — 2cr|(6). Thus, a rescattering term containing two collisions gives a negative 
contribution to the total pA cross-section. 

Let us now consider their contributions to da/dy. They are given by a\(b) + 
2a 2 (b) = 0. Indeed, in the case of a double inelastic collision, the triggered particle 
can be emitted in either of them - hence an extra factor 2. We see, in this way that 
the different contributions to da/dy of a double scattering diagram cancel. It is easy 
to see that this cancellation is valid order by order in the total number of collisions. 
This can also be seen as follows. The total inelastic cross-section for pA collision in 
the Glauber model is given by the well known expression 1 

= E *n(b) = 1 - (1 - o mel T A (b)) A . (9) 

n=l 

This expression contains a term in A 1 (single-scattering or Born term). It also 
contains contribution from multiple scattering with alternate signs (shadowing cor- 
rections). Numerically, it behaves as A a with a ~ 2/3. The single particle inclusive 
cross-section is given by 

da pA A 

——(6) oc n a n (b) = A a inel T A {b) . (10) 

<ty n=l 

the literature, the optical limit expression is often used instead of Eq. (JHJ), namely 

cr n (b) = {a ine i AT A(b)) n exp (-a me iATA(b)) jn\ . 

In this case we get <r? „,(&) = 1 — exp(— (7j ne ;ATyi(&)), which coincides with Eq. © in the large 
A limit. The same type of expressions are obtained in the eikonal model, since, in this case, the 
number of rescatterings is infinite. 
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We see that here multiple-scattering contributions cancel identically and only the 
Born term is left (impulse approximation). As a consequence of this cancellation the 
A-dependence of da/dy in pA interactions behaves as A 1 . In the case of AB collisions 
it behaves as AB and dN AB /dy = (1/ 'a AB)do AB j ' dy is proportional to A 4 / 3 , i.e. to 
the number of binary collisions - rather than to the number of participants, as one 
would naively expect (see Section 

2.3 Shadowing Corrections in the Inclusive Cross-Section 

The interest of this way of looking at the Glauber model resides in the fact that, 
as discussed in Section l2~2"l it provides strict relations between contributions to the 
total cross-section and contributions to various inelastic processes of multiparticle 
production, which make up the total cross-section via unitarity. These relations 
are the so-called Abramovsky-Gribov-Kancheli (AGK) cutting rules and have 
a general validity in RFT. The absence of shadowing corrections in the inclusive 
cross-section, Eq. (jlOj) . is called the AGK cancellation. As mentionned above, in the 
eikonal and Glauber models only the initial state is present in the vertex function 
(blob). Thus a secondary can only be produced in an interaction and the AGK 
cancellation is exact. In a general theory with a more complicated vertex function, 
the triggered particle can be produced in the blob. This gives rise to a violation of 
the AGK cancellation - which is responsible for the shadowing corrections to the 
inclusive spectra. Indeed, if the measured particle (trigger) is produced in the vertex 
function of a double scattering diagram, the extra factor 2 in c^ip) is not present 
and the AGK cancellation is not valid. In this case the shadowing corrections in 
the inclusive cross-section are the same as in the total cross-section. This is the 
physical origin of the AGK violations present in the microscopic model described in 
Sectional It is clear that if the blob has a small extension in rapidity, production 
from the blob will mainly contribute to the fragmentation region. Therefore, at 
mid-rapidities, and large energy, the AGK cancellation will be valid. 

Let us consider next the contribution to the total cross-section resulting from 
the diffractive production of large mass states. Clearly, this is equivalent to an 
increase of the rapidity extension of the blob - which, in this case, can cover the 
mid-rapidity region. Therefore, shadowing corrections to the single particle inclusive 
cross-section can be present, at mid-rapidities, provided the measured particle is part 
of the diffractively produced system. We see in this way that shadowing corrections 
to da/dy are related to diffractive production of large mass systems. The theoretical 
expression of the diffractive cross-section is well-known. It has also been measured 
experimentally and, thus, the shadowing corrections can be computed with no free 
parameters. In these lectures I will not elaborate any further on this last point. 
More details, as well as numerical calculations, can be found in |14j . 
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2.4 Space-Time Development of the Interaction : Absence 
of Intra-Nuclear Cascade 

The reggeon calculus or reggeon field theory (RFT) ^Tj provides the field theoretical 
formulation of the eikonal (for hh collisions) or the Glauber (for hA and AB) models, 
valid at high energies. The main difference between the RFT and the Glauber model 
is that, at high energies, the coherence length is large and the whole nucleus is 
involved in the interaction. Moreover, due to the space-time development of the 
interaction, when, at high energy, a projectile interacts inelastically with a nucleon 
of the nucleus, the formation time of (most of) the produced particles is larger than 
the nuclear size and, thus, most particles are produced outside the nucleus. Only 
slow particles in the lab reference frame are produced inside the nucleus and can 
interact with the nucleons of the nucleus they meet in their path (intra-nuclear 
cascade). At high energies, most of the produced particles have left the nucleus 
at the time they are formed. This near absence of nuclear cascade is well known 
experimentally. Actually, it constituted for a long time one of the main puzzles of 
high energy hadron-nuclear interactions. 

Another consequence of the space-time development of the interaction, is that 
planar multiple-scattering diagrams give a vanishing contribution to the total cross- 
section at high energies. Indeed, as discussed above, the formation time of the 
multiparticle state is larger than the nuclear size and, therefore, there is no time for 
its rescattering with other nucleons of the nucleus. The relevant multiple scattering 
diagrams are non-planar ones, describing the "parallel" interactions of different con- 
stituents of the projectile with the target nucleons (in the case of an hA collision). 
This picture is in clear contrast with the Glauber model, in which the projectile 
undergoes successive (billiard ball type of) collisions with the nucleons of the target. 

In spite of these differences, one recovers the Glauber formula in first approx- 
imation. As discussed above, this formula corresponds to the contribution of the 
initial state (on-shell projectile pole) to the various rescattering terms. In RFT one 
has, besides these contributions, also the contributions due to low mass and high 
mass diffractive excitations of the projectile. The latter are very important since, as 
we have seen in Section 12.31 they give rise to shadowing corrections to the inclusive 
cross-section. 

2.5 Regge Poles : The Pomeron 

Regge poles ^5] PZIj play an important role in high-energy interactions. It is impor- 
tant to give some basic concepts on Regge poles. Indeed, they provide a theoretical 
basis for the parametrization of the Born term in multiple-scattering models. More 
precisely, they give an important connection between high-energy behaviour and the 
spectrum of resonances in the t-channel. Also, it is very important for our purpose 
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in these lectures that they allow to determine, to a large extent, the momentum dis- 
tribution functions of partons in hadrons, as well as their fragmentation functions. 
These are the main ingredients of the microscopic models introduced in the next 
section. 

At high energies the exchange of a particle of spin J and mass Mj in the t-channel 
of a two-body amplitude 1 + 2 — > 3 + 4 has the form 



When J is large (J > 1), the cross-section increases as a power s J ~ 1 and violates 
the Froissart bound. This problem provided one of the main motivations for the 
introduction of Regge poles. Note that Eq. (jTTj) is valid only close to the pole, 
i.e. for t ~ Mj. However, the physical values of t for the s-channel amplitude 
1 + 2 — > 3 + 4 are t < 0. In the Regge pole model, the behaviour (fTT]) is modified 
for physical values of t, as follows 



Here g(t) are called Regge residues and a(t) is a function, called Regge trajectory, 
such that a{t = Mj) = J. Thus, if a(t) < 1 for physical values of t, the Froissard 
bound will be satisfied, irrespective of the value of J. For practical purposes one 
proceeds as follows. Let us consider the existing particles and/or resonances with the 
conserved quantum numbers (isospin, parity etc) of the t-channel. For each of them, 
let us plot its spin versus its mass (see Fig. 4). This is the so-called Chew-Frautschi 
plot jlbj . (Actually, in a relativistic theory it is necessary to consider separately even 
and odd values of spin). It turns out that a large family of resonances, including p, ui, 
f and A2 lie practically in a single straight line an(t) = an(0)+a't. Other resonances 
as K* and lie on parallel Regge trajectories with a(t) < au(t). Extrapolating to 
the physical region t < one finds an intercept «r(0) ~ 1/2 and a slope a' R ~ 
0.9 GeV -2 . Since a^t) < 1/2 for t < 0, the Froissard bound is respected. The 
function i](a(t)) = —[1 + <rexp(— ma(i))]/ sin7ra(t) with a = +1 (—1) for J even 
(odd), a is called the signature : positive (negative) for a = +1 (— 1). Note that for 
a trajectory of positive (negative) signature, the numerator of rj cancels the poles in 
the denominator corresponding to odd (even) values of spin. Note also that r](a(t)) 
determines the phase of the amplitude in terms of a(t). In this way the phase of the 
amplitude is related to its high energy behaviour - also determined by a(t). (The 
validity of this relationship is much more general than the Regge pole model). 

The Regge pole model has many features that have been tested by experiment 
(and no contradiction with experiment has been found). Among its successes is the 
factorization of Regge residues, Eq. (fT2"|) . Another success is that, for any given two 
body process, when resonances exist in its t-channel, a (shrinking) peak is found in 



T(s,t) = g 13 g 24 s J /(Mj-t) . 



(11) 



T(s,t) = g l3 (t) g 24 (t) s a ^ r,(a(t)) . 



(12) 
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the s-channel amplitude near t = - as obtained from Eq. (jl2J) 2 . On the contrary, 
if the t-channel is "exotic" (i.e. no known resonance exists having its quantum 
numbers), the forward peak in the s-channel is absent. There is no exception to this 
rule. 

Unfortunately, there is an important caveat. All Regge trajectories correspond- 
ing to known resonances have an intercept a(0) < 1. Therefore, it is not possible 
to explain the increase with energy observed in total cross-sections. In order to do 
so, it is necessary to postulate the existence of a Regge pole with vacuum quantum 
numbers (so that it can be exchanged in all elastic amplitudes) and intercept larger 
than one. It is called the Pomeron and is believed to correspond to the exchange of 
glue-balls in the t-channel. To compute its trajectory in QCD is a very difficult task. 
Both non-perturbative and perturbative methods give indications that its trajectory 
is slightly above one (ap(0) ~ 1.1 -v- 1.3) as required by experiment. In this way it 
is possible to explain the increase with energy of total cross-sections, i.e. a to t oc s A 
with A = a P (0) - 1 > 0. 

Note that for ctp(0) = 1, the amplitude at t = is purely imaginary. With «p(0) 
slightly above unity the ratio of real to imaginary parts is small. Unfortunately, there 
is still a caveat. Since A = ap(0) — 1 > 0, the Pomeron exchange again violates the 
Froissard bound. (This violation was one of the main motivations for Regge models, 
as discussed at the beginning of this subsection). The only way out is to use the 
Pomeron contribution as a Born term in a unitarization scheme - such as the eikonal 
or Glauber models discussed above. Indeed, if the contribution of the Born term 
or Pomeron (Fig. la) to o tot behaves as s A , the sum of the series (Figs, lb and lc) 
behaves as s A with A' < A. Moreover, in the limit s — > oo the power behaviour of 
the Born term is converted into (in s) 2 , i.e. the maximal energy growth allowed by 
the Froissard bound 3 . 

Technically, Regge poles are isolated poles in the complex angular momentum 
plane. Their s-channel iteration (as in the eikonal model) gives rise to cuts in this 
plane (Regge cuts). The latter violate factorization. However, it turns out that the 
sum of all eikonal diagrams is approximately factorizable. 

2 Assuming an exponential ^-dependence for g(t), we obtain from Eqs. and (|12|l . T(s,t) oc 
exp[(2_B + a' in s)i\. We see that the forward peak has a width that increases logarithmically 
with increasing energy. This is called the shrinking of the forward peak and has been observed 
experimentally In this way r oc in s, i.e. the effective radius of the hadron increases like in s as 
s — ► oo 

3 This can be seen as follows. If the cross-section of the Born term tends to infinity the inelastic 
cross-section of the eikonal sum, at fixed b, tends to 1 (see footnote 1). Upon integration in b one 
obtains a geometrical cross-section, proportional to r 2 oc (in s) 2 (see footnote 2). 
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3 Microscopic String Models 



3.1 Hadron-Hadron Interactions 

The Dual Parton Model (DPM) [7] |S] and the Quark Gluon String Model (QGSM) 
jU] JO] are closely related dynamical models of soft hadronic interactions. They 
are based on the large-iV expansion of non-perturbative QCD 4 and on Gribov's 
Reggeon Field Theory (RFT) [TJ]. Their main aim is to determine the mechanism of 
multiparticle production in hadronic and nuclear interactions. The basic mechanism 
is well known in e + e~ annihilation (Fig. 5). Here the e + e~ converts into a virtual 
photon, which decays into a qq pair. In the rest system of the virtual photon the 
quark (colour 3) and the antiquark (colour 3) separate from each other producing 
one string (or chain) of hadrons, i.e. two back-to-back jets. Processes of this type 
are called one-string processes. 

In hadron-hadron interactions, a one-string mechanism is also possible but only 
in some cases, namely when the projectile contains an antiquark (quark) of the same 
type than a quark (antiquark) of the target, which can annihilate with each other 
in their interaction. For instance in n + p, the d of tt + can annihilate with the d of 
p and a single string is stretched between the u of tc + (colour 3) and a diquark uu 
of p (colour 3). This mechanism is also possible in pp interactions (Fig. 6) but not 
in pp. This already indicates that it cannot give the dominant contribution at high 
energy. Indeed, when taking the square of the diagram of Fig. 6 (in the sense of 
unitarity) we obtain a planar graph, which is the dominant one according to the 
large-iV expansion. However, this only means that this graph has the strongest 
coupling. Since flavour quantum numbers are exchanged between projectile and 
target, this graph gives a contribution to the total cross-section that decreases as 
an inverse power of s (l/yfs). A decrease with s is always associated with flavor 
exchange. For instance, the charge exchange n~p — > 7r°n cross-section also decreases 
as 1/y/s. As we have discussed in Section l2*3| only an exchange in the t-channel with 
vacuum quantum numbers (Pomeron), gives a contribution to Otot which does not 
vanish asymptotically. Actually, the diagram in Fig. 6 corresponds to the exchange 
of a Reggeon, with intercept ckr(O) ~ 1/2. 

In order to prevent the exchange of flavour between projectile and target, the d 
and d have to stay, respectively, in the projectile and target hemispheres. Since they 
are coloured, they must hadronize stretching a second string of type d-d. We obtain 
in this way a two-string diagram (Figs. 7-9)). Taking the square of this diagram, 
we obtain a graph with the topology of a cylinder (Fig. 10). It turns out that this 

4 The Feynman graphs of a gauge field theory with N degrees of freedom can be classified 
according to their topology. The graphs with the simplest topology are dominant. The contribution 
of graphs with more complicated topology (characterized by well defined topological indices) are 
suppressed by powers of 1/N [17-19]. 
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is the simplest topology one can construct which does not vanish as s — > oo due to 
flavour exchange. Therefore, we obtain in this way the dominant graph for hadron- 
hadron scattering at high energy. The diagram in Fig. 10 corresponds to a Pomeron 
(P) exchange and the graphs in Figs. 7-9 are called a cut Pomeron. Its order in 
the large-iV expansion is 1/N 2 [17-19]. Note that due to energy conservation the 
longitudinal momentum fractions taken by the two components at the string ends 
have to add up to unity. 

There are also higher order diagrams (in the sense of the large- N expansion) with 
4, 6, 8 strings which give non- vanishing contributions at high energy. An example 
of the next-to-leading graph for pp interactions is shown in Fig. 9. It contains four 
strings - the two extra strings are stretched between sea quarks and antiquarks. 
The square of this graph corresponds [2U] to a two Pomeron exchange (Fig. lb) and 
has the topology of a cylinder with a handle. Its order in the large-iV expansion 
is 1/N 4 . The one with six strings corresponds [20] to a three Pomeron exchange 
(Fig. lc) and to the topology of a cylinder with two handles (order 1/N 6 ), etc. 

In this way, the large- N expansion provides the microscopic (partonic) descrip- 
tion of the Reggeon and Pomeron exchanges and of their s-channel iterations (Figs. 3), 
which were discussed in Section [2] in a hadronic picture. 

The single particle inclusive spectrum is then given by [7] 

dN pp 1 / 

-r-{y) = ^E ff n (N^(y)+N^(y) + (2n-2)N^(y) 

n 

~ iVf^(y) + Nl^{y) + {2k - 2)N^{y) (13) 
where k = J2 na n / a n is the average number of inelastic collisions. Note that each 

n n 

term consists of 2n strings, i.e. two strings per inelastic collisions. Two of these 
strings, of type qq-q, contain the diquarks of the colliding protons. All other strings 
are of type q-q. 

The weights a n of the different graphs, i.e. their contribution to the total cross- 
section, cannot be computed in the large- N expansion. However, as discussed above 
there is a one-to-one correspondence [201 between the graphs in the large- iV expan- 
sion and those in a multiple scattering model (Figs. 1). Thus, we use the weights 
obtained from the latter - with the parameters determined from a fit to total and 
elastic cross-sections (see Section |2J). At SPS energies we get k = 1.4 and at RHIC 
k = 2 at v/i = 130 GeV and k = 2.2 at y/s = 200 GeV [21]. 

The hadronic spectra of the individual strings N(y) are obtained from convo- 
lutions of momentum distribution functions, giving the probability to find a given 
constituent (valence quark, sea quark of diquark) in the projectile or in the target, 
with the corresponding fragmentation functions. Let us consider, for instance, one 
of the two qq-q strings in Fig. 9. As shown in this figure, the total energy y/s in the 
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pp center of mass frame (CM) is shared between the two strings. If \fs str denotes 
the invariant mass of a string, we have s s t r = 5x2(1 — %i) where 1 — x\ = x + and 
X2 = X- are the light-cone momentum fractions of the constituents at the string 
ends, qq and q, respectively. For massless quarks, the rapidity shift between the pp 
CM and the CM of the string is A = 1/2 log(x + /x_). We then have 

r i r i dN qq ~ q 
Nr q (s,y) = dx + dx„ pf (x + )pl(x.)— (y - A; s str ) . (14) 

Jo Jo ay 

The subscript 1 in iV and p indicates that there is only one interaction (two strings) . 
The momentum distribution functions p give the probability to find a quark or di- 
quark in the proton carrying a given momentum fraction. dN qq ~ q / dy is the rapidity 
distribution of hadron h in the CM of the qq-q string obtained from q and qq frag- 
mentation functions : 

dN qq ~ q (y - A; s str ) = f x h D qq ^ h {x h ) y > A , 

dy \ x h D q ^ h (x h ) 2/ < A , ' 

where 

1 /2 

x h = \2fi h smh(y- A)/y/s^\ , x h = (x 2 h + Ap\/ s str ^j . (16) 

Ph is the transverse mass of the detected particle h, and D q _^ h and D qq ^ h are the 
quark and diquark fragmentation functions. Momentum distribution and fragmen- 
tation functions can be obtained from Regge intercepts. Let us discuss first the 
former. In order to determine the behaviour near x = of the momentum distri- 
bution of a quark in a proton, it is convenient to look at the diagram in Fig. 6. As 
discussed in Section 12*31 the square of this diagram (in the sense of unitarity) gives a 
contribution to the total cross-section s ^ ^ 1 = e Ay ^ aR ^ where Ay = y — y ma x- 
Here y is the quark rapidity and y ma x its maximal value. Recalling that dy = dx/x, 
we obtain p q {x q ) oc x~ aR ^ = 1 — > 0. In order to determine its behaviour 
as x q — > 1, we have to use the momentum conservation x q + x qq = 1 (see Fig. 9). 
Thus, in order to have x q — > 1 it is necessary that x qq — > 0. The corresponding Regge 
exchange in the t-channel consists of two quarks and two antiquarks. Such a state is 
called a baryonium and the corresponding Regge intercept is known experimentally 
to be —1.5 ± 0.5. Taking the product of x — > and x — ► 1 behaviours we obtain 



p\{x q ) = pf{x m ) = Cx q ^ 2 xlf 5(l-x q - x qq ) = G^—{\ - x 9 ) L5 (17) 

v 1 

C is a constant determined from the normalization to unity. We see from Eq. (J17|) 
that, in average, the quark is slow and the diquark fast. 
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In order to generalize Eq. (JT7j) to the case of n inelastic interactions (2n strings), 
we just take the product of factors giving the x — > behaviour of each constituent, 
times a 5-function of momentum conservation. The momentum distribution func- 
tion p n (x) of each individual constituent is then obtained by integrating over the 
x- values of the other 2n — 1 constituents 5 . In this way the behaviour x — > is 
unchanged, whereas the power of 1 — x increases with n, due to momentum conser- 
vation. Indeed, the average momentum fraction taken by each constituent decreases 
when the number of produced strings increases. Obviously in the case of n inelastic 
collisions Eq. is still valid with p\ replaced by p n . All details can be found in 



The same Regge model considerations allow to determine the Xh — ► 1 behaviour 
of the fragmentation functions. Writing 



for the fragmentation function of constituent i into hadron h (see Eq. (jl5j)). one finds 
[22] Pi = — otfcfe(O) + A where 0^(0) is the intercept of the (kk) Regge trajectory, k 
is the system of leftover (spectator) constituents and A is a constant resulting from 
transverse momentum integrations and estimated to be A ~ | 22,. For example, 
for the fragmentation u — > 7r + , the system k is a <i-quark and j3^ + = — a^(0) + 
A. Likewise, for the fragmentation of a ud diquark into a proton we have (5 v ud = 
— a„u(0) + A and for that of a ud diquark into A, f3^ d = — a sS (0) + A. Here a U u(0) = 
a d(j(0) = a R(0) = 1/2 and a sS (0) = ot^(0) = 0. This gives a different behaviour of 
the p and A inclusive spectrum which is observed experimentally. 

Note that in writing Eq. (J13)) we have assumed that individual strings are in- 
dependent from each other. In this way, the hadronic spectra of a given graph are 
obtained by adding up the corresponding ones for the individual strings. This leads 
to a picture, in which, for any individual graph, particles are produced with only 
short-range (in rapidity) correlations. Long-range correlations (and a broadening of 
the multiplicity distributions) are due to fluctuations in the number of strings, i.e. 
to the superposition of different graphs with their corresponding weights. This gives 
a simple and successful description of the data in hadron-hadron and hadron-nucleus 
interactions [7-10]. 

3.2 Nucleus-Nucleus Interactions 

The generalization of Eq. (jSJ) to nucleus-nucleus collisions is rather straighforward. 
For simplicity let us consider the case of AA collisions and let ha and n be the 

5 Taking the same 1/yfx behaviour for both valence and sea quarks [H] JHjj these integrals can 
be performed analytically and one gets : p^{x) = C n x~ l / 2 (\ — x) n+1 / 2 and p^ q (x) = C' n x l - b (l — 
a; )n-3/2 ) with Cn = r( n + 2)/r(l/2)r(n + 3/2) and C' n = T{n + 2)/r(5/2)r(n - 1/2). 



17J 9, 




(18) 
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average number of participants of each nucleus and the average number of binary 
NN collisions, respectively 6 . At fixed impact parameter b, we have [23] [Zj 



dN AA , 
dy 



(b) = n A (b) [Nffi*(y) + N^\y) + (2k - 2)N*tf>(y) 



+(n(b)-n A (b))2kN^{y). (19) 

The physical meaning of Eq. (|19|) is quite obvious. The expression in brackets 
corresponds to a NN collision. Since n A nucleons of each nucleus participate in the 
collision, this expression has to be multiplied by n A . Note that in Eq. (fTH|) the 
average number of collisions is k - and the number of strings 2k. In the present 
case the total average number of collisions is kn - and the number of strings 2kn. 
The second term in Eq. (J2J) is precisely needed in order to have the total number 
of strings required by the model. Note that there are 2n A strings involving the 
valence quarks and diquarks of the participating nucleons. The remaining strings 
are necessarily stretched between sea quarks and antiquarks. The value of /j,(b) is 
given by //(&) = ku(b) with u(b) = n(b)/n A (b), fi(b) represents the total average 
number of inelastic collisions suffered by each nucleon. Actually, Eq. (|T9*j) is an 
approximate expression, involving the same approximation as in Eq. ([130. The 
exact expression can be found in [23] jZj- 

We see from Eq. (119)1 that dN AA /dy is obtained as a linear combination of 
the average number of participants and of binary collisions. The coefficients are 
determined within the model and depend on the impact parameter via fi(b). As 
discussed in Section l3~T1 the average invariant mass of a string containing a diquark 
at one end is larger than the one of a q-q string since the average momentum fraction 
taken by a diquark is larger than that of quark. It turns out that the same is true 
for the central plateau, i.e. : N qq ~ q (y* ~ 0) > N q ~ q (y* ~ 0). Let us now consider 
two limiting cases : 

dN AA 

If N qs ~ q °(y* ~ 0) < N qq - qv (y* ~ 0) , then — — (y* ~ 0) ~ n A ~ A 1 (20) 

dy 



If N q °- q °(y* ~ 0) ~ N qq - q "(y* ~ 0) , then — — (y* ~ 0) ~ n ~ A 4/3 . (21) 
dy 

6 n A (b) = J d 2 sAT A (s) [1 - exv{-a pp AT A {b - a)] /a AA (b) and n(b) = a pp J d 2 sA 2 T A {s)T A (b - 
s)]/a AA (b) = a pp A 2 T AA (b)/a AA (b). These expressions can be obtained in the Glauber model as 
follows. One has to generalize Eq. (JHJ to the case of AB collisions. The corresponding cross-sections 
&n A ,n B ,n{b) depend on three indices : n A (riB) is the number of participants of nucleus A(B) and 
n is the number of NN collisions. Then, n A (b) = n A a nAl n Bin (b)/ X) < J n A .n B ,n{b) and 

fry E (b). 

TtA^B »Ti TLA ,71b ,n 



16 



In the first case we obtain a proportionality in the number of participants whereas 
in the second case we obtain a proportionality in the number of binary collisions. 
Since dN AA /dy = (l/a AA )dcr AA /dy, the latter result implies that da AA /dy ~ A 2 , 
i.e. all unitarity corrections cancel and we obtain the same result as in the impulse 
approximation (Born term only). This result is the AGK cancellation discussed in 
Section|21 It implies that, for the inclusive cross-section, soft and hard processes have 
the same A-dependence. However, as discussed in Section |2] the AGK cancellation 
is violated by diagrams related to the diffraction production of large-mass states. 
These diagrams give rise to shadowing corrections. Their effect is very important in 
nuclear collisions since they are enhanced by A 1 / 3 factors. 

3.3 Charged Particle Multiplicities 

At SPS energies the limit given by Eq. (|2"Tj) is not reached, and Eq. (fHJj) leads to an 
A dependence of dN AA /dy at y* ~ in A a with a only slightly above unity, (a ~ 
1.08 between 2 and 370 participants). On the other hand, shadowing corrections 
are small due to phase space limitations. The results [21] for Pb Pb collisions at 
y/s = 17.3 GeV are shown in Fig. 10. We see that both the absolute values and the 
centrality dependence are well reproduced. When the energy increases, Eq. fl2Tj) 
shows that the value of a should increase towards 4/3, in the absence of shadowing 
corrections. However, the effect of the latter is increasingly important and, as a 
result, the value of a varies little with s. At y/s = 130 GeV, without shadowing 
corrections the A-dependence is A a , with a ~ 1.27 in the same range of n par t - a 
value which is not far from the maximal one, a = 4/3 from Eq. ()21j) . With the 
shadowing corrections taken into account, the A-dependence is much weaker (lower 
line of the shaded area in Fig. 11) |2J- The A-dependence is now A a with a ~ 1.13 
- always in the range of n part from 2 to 370. As we see, the increase of a from SPS 
to RHIC energies is rather small. This value of a is predicted to change very little 
between RHIC and LHC, where a ~ 1.1. For, the increase from a ~ 1.27 to a ~ 4/3 
obtained in the absence of shadowing is compensated by an increase in the strength 
of the shadowing corrections, leaving the effective value of a practically unchanged. 
This implies that dN/dy at y* ~ in central Au Au collisions will increase by a 
factor 2 2.5 between RHIC and LHC. This increase is slightly smaller than the 
corresponding increase of da/dy in pp collisions. 

4 Nuclear Stopping 

In pp collisions the net proton (p-p) distribution is large in the fragmentation regions 
and has a deep minimum at mid-rapidities. In contrast to this situation a much 
flatter distribution has been observed |21j in central Pb Pb collisions at CERN-SPS. 
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In view of that, several authors have claimed that the stopping in heavy ion collisions 
is anomalous, in the sense that it cannot be reproduced with the same mechanism 
(and the same values of the parameters) used to describe the pp data. In a recent 
paper [23] it has been shown that this claim is not correct. 

In the model described in the previous section, the net baryon can be produced 
directly from the fragmentation of the diquark. Another possibility is that the 
diquark splits producing a leading meson in the first string break-up and the net 
baryon is produced in a further break-up. Clearly, in the first case, the net baryon 
distribution will be more concentrated in the fragmentation region than in the second 
case. The corresponding rapidity distributions are related to the intercepts of the 
relevant Regge trajectories, a qq and a g , respectively, i.e. they are given by e A2y ( 1_a ) 
(see Section ETTj) . Here Ay is the difference between the rapidity of the produced 
net baryon and the maximal one. In the case of the first component, in order to 
slow down the net baryon it is necessary to slow down a diquark. The corresponding 
Regge trajectory is called baryonium and its intercept is known experimentally to 
be a qq = a qqqq (0) = —1.5 ± 0.5 (see Section ETT]) . For the second component, where 
a valence quark is slowed down, we take a q = ct qq (0) = «r(0) = 1/2 7 . 

In this way we arrive to the following two component model for net baryon 
production Bi — Bi (where i denotes the baryon species) out of a single nucleon 



A fjBi—Bi 

— (y,b) = P 2 a Cu(b) Z\- a ^\l - 2" )M(6)-3/2+n.,(a«(0)-^(0)) 

dy ' 

- a)C' Kb) Z+" Qw(0) X (1 - Z+ )Mb)-3/2+c+n sg (a ii (0)-a <A (0)) ^ 

where n sq is the number of strange quarks in the hyperon, «_r(0) = 1/2, a^(0) = 
0, Z + = {t v ~ Vmax \ y ma x is the maximal value of the baryon rapidity and nip) is 
the average number of inelastic collisions suffered by the nucleon at fixed impact 
parameter b (see Section EOj) . The constants C M and C are normalization constants 
required by baryon number conservation 8 . The small Z behaviour is controlled by 
the corresponding intercept. The factor (1 — Z^.)^ -3 / 2 gives the Z — > 1 behaviour 
of the diquark momentum distribution function in the case of [i inelastic collisions 
(see footnote 5). Following conventional Regge rules an extra an(0) — 0^(0) = 1/2 
is added to the power of 1 — Z + for each strange quark in the hyperon (see Section 
EU Eq. (HHD). 

7 There is a third possibility in which the net-baryon transfer in rapidity takes place without 
valence quarks (string junction or gluonic mechanism) with intercept either asj = 1/2 [2(ij or 
as j = 1 |27| . We find no evidence for such a component from the existing pp and AA data. Its 
smallness could be related to the fact that it produces an extra string of hadrons and, thus, does 
not correspond to the dominant topology in the large N expansion. 

8 CV = T(a + b)/T{a)T{b) with a = 1 - or 9 (0) and b = fi(b) - 1/2 + n sq (a R (0) - o^(0)) ; 
C' v = Y{a! + b')/T(a')T(b') with a' = 1 - a qq (0) and b' = fi(b) - 1/2 + c + n sq (a R (0) - a (O)). 
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The weights I\ and I\ allow to determine the relative yields of the different 
baryon and antibaryon species. They are computed in Appendix A using simple 
quark counting rules. 

The fraction, a, of the first component is treated as a free parameter. The 
same for the parameter c in the second component - which has to be determined 
from the shape of the (non-diffractive) proton inclusive cross-section in the baryon 
fragmentation region. It can be seen from Eq. (J22)) that stopping increases with 
//(&), i.e. with the total number of inelastic collisions suffered by each nucleon. This 
effect is present in the two terms of and is a consequence of energy conservation. 
The question is whether this "normal" stopping is sufficient to reproduce the data. 
In other words whether the data can be described with a universal value of a, i.e. 
independent of \x and the same for all reactions. 

The formulae to compute net baryon production in pp or AA collisions can be 
obtained from Eq. (J2~2"j) in a straightforward way. Thus, in AA collisions we have : 
dN AA ^ B ^/dy(y*, b) = n A (b)[dN B >-^/dy(y*, b) + dN B ^/dy(-y\ b)\. 

Note that in the formalism above, baryon quantum number is exactly conserved. 
Note also that shadowing corrections are not present. Indeed, as explained in Sec- 
tion |2J these corrections affect only the term proportional to the number of binary 
collisions, which is not present for net baryon production. 

A good description of the data on the rapidity distribution of pp — > p — p + X 
both at y/s = 17.2 GeV and y/s = 27.4 GeV is obtained from Eq. with a = 0.4, 
c=l, a q = 1/2 and a qq — — 1. The results are shown in Tabled] at three different 
energies, and compared with the data. As we see the agreement is reasonable. As 
it is well known, a pronounced minimum is present at y* = 0. There is also a 
substantial decrease of the mid-rapidity yields with increasing energy. Also, the 
mid-rapidity distributions get flatter with increasing energy since the net proton 
peaks are shifted towards the fragmentation regions. 

It is now possible to compute the corresponding net baryon production in heavy 
ion collisions and to check whether the data can be described with Eq. (f2"2"J) using 
the same set of parameters as in pp. The results for net protons {p — p) in central Pb 
Pb collisions at y/s = 17.2 GeV and central Au Au collisions at y/s = 200 GeV are 
shown in Fig. 14. We see a dramatic change in the shape of the rapidity distribution 
between the two energies, which is reasonably described by the model. Therefore, 
we conclude that there is no need for a new mechanism in AA collisions. 

5 Hyperon and Antihyperon Production 

Strange particle production, in particular, of multistrange hyperons, has been pro- 
posed as a signal of Quark Gluon Plasma formation. Flavor equilibration is very 
efficient in a plasma due to large gluon densities and low thresholds [4j. Moreover, 
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the increase of the relative yields of strange particles in central AA collisions as 
compared to pp can be understood as a consequence of the necessity of using the 
canonical ensemble in small size systems (pp) - rather than the grand canonical one. 
The exact conservation of quantum numbers in the former leads to a reduction of 
ss pair production, as compared to the latter [3] . 

An analysis of the results at SPS in the framework of the present model has been 
presented in |2*5j . In the following we concentrate on RHIC results. 

A general result in DPM is that the ratios B/h~ and B/h~ of baryon and 
antibaryon yields over negatives decrease with increasing centralities. This is easy 
to see from Eq. (jl9j) . The production from q s -q s strings scales with the number of 
binary collisions. These strings have a smaller (average) invariant mass than the qq-q 
strings and, thus, are more affected by the thresholds needed for BB pair production. 
As a consequence, the centrality dependence of B and B production will be smaller 
than the one of negatives. The effect is rather small at RHIC energies. However, it 
is sizable and increases with the mass of the produced baryon. In contrast with this 
situation, the data for A's show no such decrease and an increase is present for H 
production. Data on Q production are not yet available. However, SPS data clearly 
show a hierachy in the sense that the enhancement of baryon production increase 
with the mass (or strange quark content) of the produced baryon. 

The only way out we have found is to give up the assumption of string inde- 
pendence. Until now we have assumed that particles produced in different strings 
are independent from each other. In the following we allow for some final state 
interactions between comoving hadrons or partons. We proceed as follows. 

The hadronic densities obtained in Section 2 are used as initial conditions in 
the gain and loss differential equations which govern final state interactions. In 
the conventional derivation [2H] of these equations, one uses cylindrical space-time 
variables and assumes boost invariance. Furthermore, one assumes that the dilution 
in time of the densities is only due to longitudinal motion, which leads to a r _1 
dependence on the longitudinal proper time r. These equations can be written 

m 

r ~T~ = a ke Pk Pi - a ik Pi Pk ■ (23) 
aT kl k 

The first term in the r.h.s. of ()23j) describes the production (gain) of particles of 
type % resulting from the interaction of particles k and i. The second term describes 
the loss of particles of type i due to its interactions with particles of type k. In Eq. 
(|2l?|) pi = dNi/dyd 2 s(y,b) are the particles yields per unit rapidity and per unit of 
transverse area, at fixed impact parameter. They can be obtained from the rapidity 
densities, Eq. (HHJ), using the geometry, i.e. the s-dependence of tia and n. The 
procedure is explained in detail in jHOj where the pion fragmentation functions are 
also given. Those of kaons and baryons can be found in j3T] . These fragmentation 
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functions are obtained using the procedure sketched at the end of Section 13.11 (see 
Eq. ()18|)). a^i are the corresponding cross-sections averaged over the momentum 
distribution of the colliding particles. 

Equations (J2l?j) have to be integrated from initial time To to freeze-out time 77. 
They are invariant under the change r — > ct and, thus, the result depends only 
on the ratio Tf/r . We use the inverse proportionality between proper time and 
densities and put Tj/r = (dN/dyd 2 s(b))/pf. Here the numerator is given by the 
DPM particles densities. We take p/ = [3/irRp}(dN~ /dy) y *^ = 2 fm~ 2 , which 
corresponds to the density of charged and neutrals per unit rapidity in a pp collision 
at ,/s = 130 GeV. This density is about 70 % larger than at SPS energies. Since 
the corresponding increase in the AA density is comparable, the average duration 
time of the interaction will be approximately the same at CERN SPS and RHIC - 
about 5 to 7 fm. 

Next, we specify the channels that have been taken into account in our calcula- 
tions. They are 

nN <± fCA(E) , ttA(S) KE , ttS ^ KVL . (24) 
We have also taken into account the strangeness exchange reactions 

ttA(S) ^ KN , ttS ^ KA(E) , vrfi ^ KE (25) 

as well as the channels corresponding to (}2*4"|) and (|2*3j) for antiparticles. We have 
taken = a = 0.2 mb, i.e. a single value for all reactions in (J24"|) and (}25|) - the 
same value used in ref. |2E] to describe the CERN SPS data. 

Before discussing the numerical results and the comparison with experiment let 
us examine the qualitative effects of comovers interaction. As explained in the 
beginning of this Section, without final state interactions all ratios K/h~, B/h~ 
and B/h~ decrease with increasing centrality. The final state interactions (J2U), (|25|) 
lead to a gain of strange particle yields. The reason for this is the following. In the 
first direct reaction (j21J) we have p T > px, Pn > Pa, P-kPn PkPa- The same is 
true for all direct reaction (0J). In view of that, the effect of the inverse reactions 
f!24|) is small. On the contrary, in all reactions (|23j) . the product of densities in the 
initial and final state are comparable and the direct and inverse reactions tend to 
compensate with each other. Baryons with the largest strange quark content, which 
find themselves at the end of the chain of direct reactions ()24)1 and have the smallest 
yield before final state interaction, have the largest enhancement. Moreover, the gain 
in the yield of strange baryons is larger than the one of antibaryons since Pb > Pb- 
Furthermore, the enhancement of all baryon species increases with centrality, since 
the gain, resulting from the first term in Eq. ()23|) . contains a product of densities 
and thus, increases quadratically with increasing centrality. 
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In Fig. 15a-15d we show the rapidity densities of B, B and B — B versus h~ = 
dN~ /drj = {l/l.n)dN~ / dy |3T] and compare them with available data [32-34]. We 
would like to stress that the results for H and H were given [HI] before the data jS] . 
This is an important success of our approach. 

In first approximation, the yields of p, p, A and A yields over h~ are independent 
of centrality. Quantitatively, there is a slight decrease with centrality of p/h~ and 
p/h~ ratios, a slight increase of A/h~ and A/h~ and a much larger increase for 

5 (H)//t~ and Q {Vt)/h~. This is better seen in Figs. 16a and 16b where we plot 
the yields of B and B per participant normalized to the same ratio for peripheral 
collisions versus n part . The enhancement of B and B increases with the number 
of strange quarks in the baryon. This increase is comparable to the one found at 
SPS between pA and central Pb Pb collisions. (In the statistical approach the 
enhancement of B and B relative to pp decreases with increasing energy. This may 
allow to distinguish between the two approaches). 

The ratio K~/n~ increases by 30 % in the same centrality range, between 0.11 
and 0.14 in agreement with present data. The ratios B/B have a mild decrease with 
centrality of about 15 % for all baryon species - which is also seen in the data. Our 
values for N ch /N± x = 1/2 are : p/p = 0.69, A/A = 0.74, S/S = 0.79, Q/tt = .83, 
to be compared with the measured values j3SJ : 

p/p = 0.63 ± 0.02 ± 0.06 , A/A = 0.73 ± 0.03 , S/3 = 0.83 ± 0.03 ± 0.05 . 

The ratio K + /K~ = 1.1 and has a mild increase with centrality, a feature also seen 
in the data. 

Note that a single parameter has been adjusted in order to determine the absolute 
yields of BB pair production, namely the p one - which has been adjusted to the 
experimental p value for peripheral collisions. The yields of all other BB pairs has 
been determined using the quark counting rules given in Appendix A. 

Although the inverse slopes ( "temperature" ) have not been discussed here, let us 
note that in DPM they are approximately the same for all baryons and antibaryons 
both before and after final state interaction - the effect of final state interaction on 
these slopes being rather small 

6 J/ip Suppression 

A most interesting signature of the production of QGP is the suppression of res- 
onance production As a consequence of deconfinement, the resonances are 
"melted", i.e. the bound state cannot be formed. More precisely, as in the case 
of an ordinary plasma, the potential Vo(r) is screened (Debye screening), i.e. it is 
changed into V (r) exp(— r/r D (T)). Here r D (T) is the Debye radius that decreases 
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with increasing temperature. When 7\d(T) becomes smaller than the hadronic ra- 
dius, the bound state cannot be formed. This idea is particulary interesting in the 
case of the J ftp (a resonance consisting of a charm quark and its antiquark). Indeed, 
it has been shown that the melting of the J /if) occurs at temperature only slightly 
higher than the critical temperature at which QGP is formed. Moreover, the pro- 
duction of c-c pairs is very rare. If they cannot bind together the c(c) will combine 
with a light quark (antiquark) giving rise to a D(D) meson (open charm) 9 . 

The NA38-NA50 collaborations have observed a decrease of the ratio of J /if) to 
dimuon (DY) cross-sections with increasing centrality in SU and Pb Pb collisions 
|3*7] . The same phenomenon has been observed in pA collisions with increasing values 
of A. In this case, it is interpreted as due to the interaction of the pre-resonant cc pair 
with the nucleons of the nucleus it meets in its path (nuclear absorption). [Indeed, 
the formation time of the J / if) is longer and it is produced outside the nucleus.] As a 
result of this interaction, the cc pair is modified in such a way that, after interaction, 
it has no projection into J /if) (a DD pair is produced instead). The corresponding 
cross-section is denoted <r a b s (absorptive cross-section). 

The survival probability, S a b s (b), of the J /if) in pA collisions can be easily calcu- 
lated in the probabilistic Glauber model. One has 

r+oo f r+oo \ 

S±(b) = J dZ PA (b, Z) (l - a abs j p A (b, Z 1 )dZ 1 ) 

J ^ dZp A (b, Z) = {1 " [1 " °a bs T A (b)]} . (26) 

Indeed, the cc pair is produced at a point of coordinates (b, Z) inside the nucleus, 
with probability proportional to p A (b,Z). The term inside the parenthesis gives 
the probability of non-absorption during its subsequent propagation through the 
nucleus. Note that S a b s = 1 for a a b s = 0. The generalization of (J26|) to the case of 
nucleus-nucleus interactions is rather straightforward. We have [3U] 

S^(b,s) = S* bs (s) SZs(b~s) . (27) 

The NA50 collaboration has shown that the J / if) suppression in Pb Pb collisions 
has an anomalous component, i.e. it cannot be reproduced using nuclear absorption 
alone [SZj. Two main interpretations have been proposed : deconfinement and 
comovers interaction. The latter mechanism has been described in Section 03 for 
hyperon production. In the case of J /if) suppression, a single channel is important 
namely cc (or J /if)) interacting with comoving hadrons and producing a DD pair. 
In this case, Eq. (|2^jl can be solved analytically. One obtains for the expression of 
the survival probability S co |3*U] 

9 In what follows we disregard the possibility of c and c recombination into a J/ip. 
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S* B (b,s) = exp 



-<7coW DT (b,s)£n^ 

Z Tq 



(28) 



Here Ny° DT (b, s) is the density of charged particles in the rapidity region of the 
dimuon trigger (DT) (0 < y* < 1) computed from Eq. The factor 3/2 takes 

care of the neutrals. For tj/tq in Eq. (|25j) we use the expression given in Section 
with = 1.15 fm~ 2 at SPS energies. Note that S co depends on a single parameter 
<j co , the effective cross-section for comovers interaction. 

The results [38 j of the comovers interaction model are presented in Fig. 17. The 
agreement with the data |3Z| is quite satisfactory. There is a single free parameter 
a co = 0.65 mb. The value of a a b s = 4.5 mb is determined jSH] from the pA data and 
the absolute normalization (47) from the SU ones. 

Predictions of the comovers model [3B] at RHIC energies are given in Figure 18. 

In a deconfming approach one proceeds as follows [HH]- One assumes that the 
energy density of the produced system is proportional to the density of participants 
n,A(b, s) of nucleus A in a AA interaction. If n^ib, s) < n crit the J /if) is suppressed 
only due to ordinary nuclear absorption with cross-section a a b s . On the contrary, if 
n A(b, s) > ncrit, the nuclear absorption formula is used with a a b s infinity. In this way, 
no J/ if) can survive above the critical density. The deconfming approach leads to a 
satisfactory description of J /if) suppression, with n crit treated as a free parameter 10 . 
However, a quantitative analysis of the most recent NA50 data jH7j is still missing. 
On the other hand, the centrality dependence of the average px of J /if) is better 
described in the comovers approach than in a deconfming scenario 



7 Conclusions 

In these lectures a description of microscopic string models of hadronic and nucleus 
interactions has been presented. Consequences of the model for charged particles 
multiplicities, net baryon production (stopping), hyperon and antihyperon produc- 
tion (strangeness enhancement) and J /if) suppression have been examined. 

As a starting point we have assumed that particles produced in different strings 
are independent, i.e. there is no "cross-talk" between strings. While this assumption 
works quite well in pp and pA interactions, the AA data can only be described if final 
state interaction (comovers interaction) is introduced. However, the corresponding 
cross-section turns out to be rather small (a few tenths of a milibarn). Due to this 
smallness and to the short duration time of final state interaction (5 4- 7 fm) it is 
unlikely that thermal equilibrium can be reached. 

10 Both in the comovers model and in the deconfming one, the description of the J/ip suppression 
at very large transverse energy, Et, requires the introduction of £/r-fluctuations . 
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Of course it is not possible to reach a definite conclusion on this important point. 
However, particle abundances not only do not allow to conclude that equilibrium 
has been reached, but, on the contrary, their centrality dependence tends to indicate 
that this is not the case. Let us consider for instance p and p production. In our 
model, their yields are practically not affected by final state interaction, i.e. they 
are practically the same assuming string independence. Yet, the model reproduces 
the data, from very peripheral to very central interaction. This success would be 
difficult to understand in a QGP scenario in which for peripheral collisions (below 
the critical density) there is strong, non-equilibrated, pp annihilation, which becomes 
equilibrated for central ones, above the critical density. More generally, the QGP 
scenario would be strongly supported if some kind of threshold would be found in the 
strange baryon yields around the critical energy density. At SPS energies, evidence 
for such a threshold in the S yield was claimed by the NA57 collaboration based on 
preliminary data jH], but it is not seen in the more recent analysis 02]- Moreover, 
the saturation at large centralities of the B and B yields per participants, shown 
by the WA97 data, has also disappeared from the new data [42J, in agreement with 
the predictions of the comovers interaction model (SHI- Unfortunately, these data 
only cover a limited range of centrality. In contrast to this situation, the RHIC data 
explore the whole centrality range from very peripheral to very central collisions and 
the centrality dependence of the yields of p, A, 5 and their antiparticles shows no 
structure whatsoever. Moreover, the yields of H and H per participant (as well as 
the ratios H//i~ and H//t~) do not seem to saturate at large centralities. If the same 
happens for Q and Q production (as predicted in our approach, Figs. 16) the case 
for QGP formation from strange baryon enhancement will be quite weak. 

Finally, it should be stressed that the final state interaction of comovers in our 
approach is by no means a trivial hadronic effect. Indeed, the interaction of comovers 
starts at the early times when densities, as computed in DPM, are very large. In 
this situation the comovers are not hadrons (there are several of them in the volume 
normally occupied by one hadron, and, moreover, at these early times hadrons are 
not yet formed). This is probably the reason why in our approach the comover 
interaction cross-sections required to describe the data are smaller than in hadron 
gas models where the final state interaction starts only after hadron formation. 
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Table 1: Calculated values [23] of the rapidity distribution of pp — > p — p + X at 
y/s = 17.2 GeV and 27.4 GeV (k = 1.4) and y/s = 130 GeV (k = 2). (In order to 
convert da/dy into dN/dy a value of a = 30 mb has been used). For comparison with 
the nucleus- nucleus results, all values in this table have been scaled by = 175. 



y* pp —> p — p pp — > p — p pp — ► p — p 
y/j = 17.2 GeV = 27.4 GeV = 130 GeV 






9.2 


6.5 


3.6 






[6.3 ±0.9] 




1 


15.0 


9.3 


4.2 




[16.1 ±1.8] 


[9.6 ±0.9] 




1.5 


25.8 


14.6 


5.1 




[24.1 ± 1.4] 


[15.4 ±0.9] 




2 


47.1 


26.2 


6.8 




[45.4 ± 1.4] 


[27.7 ±0.9] 





Appendix A 

In order to get the relative densities of each baryon and antibaryon species we use 
simple quark counting rules [2E| [3*Tj . Denoting the strangeness suppression factor 
by S/L (with 2L + S = 1), baryons produced out of three sea quarks (which is the 
case for pair production) are given the relative weights 

J 3 = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 (A.l) 

for p, n, A±E, 2°, H~ and Q, respectively. The various coefficients of J3 are obtained 
from the power expansion of (2L + S) 3 . 

For net baryon production, we have seen in Section0]that the baryon can contain 
either one or two sea quarks. The first case corresponds to direct diquark fragmen- 
tation described by the second term of Eq. The second case corresponds to 
diquark splitting, described by the first term of ([22)1 . In these two cases, the relative 
densities of each baryon species are respectively given by 

h = L:L:S (A.2) 

for p, n and A ± E, and 

I 2 = 2L 2 : 2L 2 : ALS : ^S 2 : ^S 2 (A.3) 

for p, n, A± E, H° and The various coefficients in (A.2) and (A.3) are obtained 
from the power expansion of (2L + S) and (2L + S) 2 , respectively. 
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In order to take into account the decay of S*(1385) into An, we redefine the 
relative rate of A's and S's using the empirical rule A = 0.6(S + + £~) - keeping, 
of course, the total yield of A's plus S's unchanged. In this way the normalization 
constants of all baryon species are determined from one of them. This constant, 
together with the relative normalization of K and n, are determined from the data 
for very peripheral collisions. In the calculations we use S = 0.1 (S/L = 0.22). 
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Figure captions 

Figure 1 : Phase diagram in statistical QCD pQ 

Figure 2 : Energy density versus temperature in lattice QCD pQ at // = 0, for 
three light flavors (upper), two light and one heavier (middle) and two light flavors 
(lower) . 

Figure 3 : Single (a), double (b), and multiple (c) scattering diagrams in the eikonal 
model. 

Figure 4 : The secondary Regge trajectories with highest intercept. 

Figure 5 : The mechanism of particle production in e + e~ annihilation. The net of 
soft gluons and quark loops is only shown here and in Fig. 10. 

Figure 6 : One string diagram in pp. 

Figure 7 : Dominant two-chain (single cut Pomeron) contributions to high energy 
7r + -proton collisions. 

Figure 8 : Dominant two-chain contribution to proton-antiproton collisions at high 
energies (single cut Pomeron). 

Figure 9 : Dominant two-chain diagram describing multiparticle production in high 
energy proton-proton collisions (single cut Pomeron). 

Figure 10 : Single Pomeron exchange and its underlying cylindrical topology. This 
is the dominant contribution to proton-proton elastic scattering at high energies. 

Figure 11 : Two cut Pomeron (four-chain) diagram for proton-proton collisions. 

Figure 12 : The values of dN ch /dy per participant for Pb Pb collisions at \fs = 
17.3 GeV computed [21] from Eq. ([T9"]). compared with WA98 data. 

Figure 13 : The values of dN ch /dr] cm /(0.5 n part ) for Au Au collisions at y/s = 
130 GeV computed [21] from Eq. ()19j) including shadowing corrections are given 
by the dark band in between solid lines. The PHENIX data are also shown (black 
circles and shaded area). 

Figure 14 : Rapidity distribution of net protons (p-p) for the 5 % most central Pb 
Pb collisions at SPS (-^i = 17.2 GeV) and for the 10 % most central Au Au ones 
at ^ = 200 GeV, compared to data [H] 

Figure 15 : (a) Calculated values [21] of dN/dy of p (solid line) p (dashed line), 
and p — p (dotted line) at mid rapidities, \y*\ < 0.35, are plotted as a function of 
dNh-fdrj, and compared with PHENIX data [22] ; (b) same for A and A compared 
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with preliminary STAR data [33] ; (c) same for H and H compared to preliminary 
STAR data 34 j ; (d) same for Q and fl. 

Figure 16 : Calculated values [HI] of the ratios B/n part (a) and B/n part (b), nor- 
malized to the same ratio for peripheral collisions (n part = 18), plotted as a function 
of n part . 

Figure 17 : The ratio of J / ip over DY cross-sections in Pb Pb collisions a 158 GeV/ c 
versus Et obtained [HE] in the comovers interaction model with a a b. s = 4.5 mb and 
a co = 0.65 mb. The absolute normalization is 47. The preliminary data are from 

Figure 18 : dN/dy of J/ip (times branching ratio) in Au Au collisions at yfs = 
200 GeV per nucleon and y* ~ 0, scaled by the number of binary collisions < nb >, 
versus the number of participants. The curves are obtained in the comovers model 
|3*H] with a a bs = 0, a co = 0.65 mb (upper) and a a bs = 4.5 mb, a co = 0.65 mb (lower). 
The absolute normalization is arbitrary. An extra 20 % suppression between pp and 
central Au Au is expected due to shadowing. 
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Figure 12 
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Figure 16a 
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